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ABSTRACT 

The plane strain problem of a multi-layered composite with 
parallel cracks in considered. The main objective of this paper 
is to study the interaction between parallel and collinear cracks. 
The problem is formulated in terms of a set of simultaneous 
singular integral equations which are solved numerically. The 
effect of material properties on the interaction between cracks 
is also demonstrated. 
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INTRODUCTION 


Welded and bonded structures have been observed to contain 
multiple cracks. The study of the interaction between such 
cracks has also been of considerable interest to reactor design- 
ers. The problem of a multi-layered composite containing a 
single crack was studied by Erdogan and Gupta [1-2]. The inter- 
action between multiple cracks in an isotropic medium and collin- 
ear cracks in a layered composite has been considered by Ratwani 
[3-4]. 

In the present study, the analytical methods of [1-4] have 
been extended to treat the layered composite containing parallel 
and collinear cracks. In particular, the plane strain problem 
of an elastic layer bonded to two dissimilar half-planes is 
considered. The layer medium contains one or two symmetri cal 1 y 
placed collinear flaws and one of the half-planes is assumed to 
have a single parallel flaw. The procedure, of course, can 
easily treat any composite containing n elastic layers and cracks 
located along m parallel planes. For the sake of simplicity, 
only the symmetric problem is studied here. The anti -symmetri c 
loading case can be handled in an analogous manner. 

Stress intensity factors at all the crack tips are computed. 
Their variation with respect to the crack locations, geometry 
and material of the composite are presented graphically. 
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FORMULATION OF THE PROBLEM 


Consider the plane problem, shown in Figure 1, containing 
one or two collinear cracks in each plane. The cracks are 
assumed to be located symmetri cal ly with respect to the y-axis. 
In this paper, our primary interest is in the disturbed stress 
state caused by the cracks. Hence, assuming that the overall 
stress distribution a . .° in the imperfection-free medium is 

* vJ 

known, the stress state a..^ in the cracked medium may be ex- 

* J 

pressed as 


a • = a.. 0 + a . • 

ij ij iJ 


( 1 ) 


where a., is the disturbed stress state obtained by subjecting 
the crack surfaces to the following tractions. 

- 1 ( x , h ) = - cr yy °(x,h) = p 2 ( x ) 

, C < I x | < D , 


yy 


a xy 1 ( x 5 h ) = - a°(x,h) = Pi 


xy 


xy 


a yy 3 ( x ,_h i ) = " a yy° (x, " h i ) = ^2 (x) 

9 

a xy 3 (x,-h 1 ) = - a X y° ( x > -h-j ) = q -j ( x ) 

p^(x) and q^x) satisfy a Holder condition 
ranges . 


( 2 ) 


A< | x | <B . 

in their respective 


The integral transform technique, described in detail for 
a single crack [1] and for multiple cracks [3], is used here to 
formulate the problem in terms of four unknown functions defined 
by 
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f -J ( X ) 


3 / 

3x U 2 


u 1 ) , 


f 2 (x) = 


3 f 

3x ^ 1 * * * V 2 


V 1 ) » 


( C< | x | <D , y = h ) , 


{ 3 } 

9 1 (x) = 37 (u 3" u 4 )s 9 2 (x) = v 3 _ v 4 ) ’ ( A< I X i <B , y = -h } ). 

Note that the crack surfaces are the singular surfaces across 
which the displacement vector suffers a discontinuity and the 
unknown functions define the derivatives of the crack opening 
displacements. For the sake of simplicity, the central plane 
of the elastic layer is assumed to have the cracks. The case 
when the crack lies at the interface between the bielastic media 
has been treated in detail in [2]. 


Following the procedure of [1] and [3], a set of simultane- 
ous singular integral equations of the first kind is derived, 
expressed as follows: 

! I j, S k (T)H 1 k (x.T)dT 

L i j I L 2 K ~ I 

1 +K, 

= ~ 2yTj" u P-j( x ) » C < | x | < D , i = 1 , 2 , 

( 4 ) 

2 2 6. 

/ I f -j ( t ) L ( x , t ) d t + / l gJxH-^ + M . . ( x , x ) ] d x 
L 1 j=l J L 2 k=l K T “ x 1k 

1 +K ? 

= ~2]I^ 1r< ?i ( x ) » A< | x | <B , i =1 , 2 , 

where L 1 = ( C < J x | < D ) and L 2 = ( A< | x | < B ) , and = 3-4v i for 
plane strain and = (3-v. )/(l+v. ) for generalized plane 
stress. and v. are the shear modulii and the Poisson's 

ratios, for i = 1 , 2 , denoting the elastic layer by the subscript 
1 and the half-planes by 2. The functions K . . , H-, and M • . 

•J 1J IK IK 
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are Fredholm kernels and are bounded in their respective closed 
intervals. The expressions for these Fredholm kernels are given 
as follows: 


00 s, (a) - 4 ah 

K ll (x ’ t) = l D^Ta] 


- 2 ah . / . \ , 

e sina(t-x)da 


K-, 2 ( x , t ) = K 91 ( x s t ) = 0 


21 


00 s, (a) + 4 ah . 

K 22 (x,t) = / D^ToO e”^ an s i na ( t-x ) da 


P 2 l+<- 00 


H 11^ X,T ^ ^2" y l 2X 3 o ^ X 3 s 3^ a ^ + d-Zah^s^a)]. 


-a( h+h-j ) 

D ^j~( a ) sina(T-x)da 


y 2 1 +k 


H 12< x > t ) ■ v pfp l [ X 3 S 3<“> - (' + 2“h,)s 2 (a)]- 


■a( h+h 1 ) 


~D^W) 

y 9 1 +K, <=° 


cosa( t-x ) da 


‘21 


(x ’ T) = 2X7 1 l t A 3 s 5 (a) - 0 - 2c,h l> s 4 ( “ )]- 


•a( h + h 1 ) 

Do(a) 


cosa( x-x)da 


Up 1 + k , 00 

H 22 ( x , r) - - 2A f ^ X 3 ^ 5 ^ ^ + (1 + 2ah-|)s^(a)] 


( 5 ) 


y 2 " y l ^3 0 


e -“< h+ h 1 ) 

D 2 (a) s l na( t-x ) da 


(Po-Pl ) 

L ( x , t ) = i — H • • ( x , t ) , i , j = 1,2 

1J y 2 2 (i +K i ' 
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“ 2 j [ X 3 + S 8^ a ^ + D 1 (a) D 2 (a) { X 3 s 7 ^ a ^ 


+ [Sy ( a) Sg ( a) - 8ah(l - 2ah-| )]e" 4ah } 


!~ 2ah l sina(x-x)d 


a 


M 12 (x,t) = - \ f 


( X i - X 2 ) 


3 ~ s 6 ^ a ^ + D 1 (a)D 2 (a) {^ 3 S 7 ( a ) 

e 2a ^l cosa(x-x)da 


+ [ 16 a 2 hh 1 - s 6 (a)s 7 (a}]e" 4ah } 


M 2 i (x,x) = M 12 (x,t) 
M 22 (x,t) 


i 00 r ( x i _ x 2 ) 

- 2 / [x 3 + s g (a) + (a)D 2 (a) { x 3 s 7 ( ot ) 


- [s 7 (a)s g (a) - 8ah ( 1 + 2ah-| ) ]e~ 4a ^ } 


-2ah] 


sina(x-x)da 


where 

D-j ( a) = - x 2 + ( 4 ah + X ] e" 2ah ) e _2ah 
D 2 ( a ) = X 2 + ( 4 ah- X ] e" 2ah )e' 2ah 
S-| ( a) = 1 + X^ 2 + 4 a 2 h 2 - 2 X ie " 2ah 
s 2 (a) = - X 2 + (1 + 2ah)e" 2atl 
S3 ( a) = 1 - 2ah - X ie " 2ah 
S4 (a) = X 2 - (1 - 2ah)e' 2ah 
s 5 (a) = - 1 - 2ah + X^ -2 ®* 1 
s 6 (a) = (1 - 4 a 2 h 1 2 )/X 3 
s 7 (a) = - X 2 + A-j e 


( 6 ) 
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Sg(a) = (1 - 2ah 1 ) 2 /A 3 

s g (a) = (1 + 2ah ] ) 2 /A 3 
and 


A 


1 


y-j <2 ~ y 2 K l 

u 2 + i j i k 2 


A 


2 


y 1 y 2 ^ i 
^ ' y 2 


A 


3 


y 2 V* ■] k 2 

y 2 “ 


(7) 


The unknown functions f. and in equations (4) have 
integrable singularities at the end points. Therefore, the 
equations (4) must be solved subject to the singlevaluedness 
conditions 

D B 

/ f,(t)dt = 0 = / g • ( t ) d t , 1 = 1,2. (8) 

C 1 A 1 

The singular integral equations (4) are solved simultane- 
ously by using the numerical technique described in [5]. It may 
be noted that if one of the cracks lies on the interface, the 
corresponding integral equation would become that of second kind. 
The numerical technique to treat such equations is described in 
[6] and is used to solve the interface crack problems in [2]. 

The stress intensity factors Kj and Kjj at all the crack tips 
are defined as in [1]. As an example, for the crack in medium 
(1) near the crack tip x-*D, these can be expressed as 
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( 9 ) 


] T = lim /2(x-D) o '(x.h) 
1 v+n y y 


x+D 

K^t = lim /2(x-D) o ^(x.h) . 
x+D xy 

The stress intensity factors can also be expressed in terms of 
the unknown functions f i (x) . The functions f^x), which have 
integrable singularities, may be written as 

G,(x) 

f ,• ( x ) = — - — — • 

/(D-x) ( x-C) 

Equations (9) can now be written as 

I 

K I = - t+t— lim /2 ( D-x ) f~(x) 

1 1 K 1 x+D d 


( 10 ) 


2y 


1+K 


K n /d?c 


2y 


(ID 


II = TTi h li * ^xr ^(x) 


1 x->D 


2y 


1+K 


K i /Zl G 1 


Superscripts 1 and 2 on the stress intensity factors refer to 
the cracks in medium 1 and 2 respectively. 


DISCUSSION OF RESULTS 

To demonstrate the interaction between parallel cracks, a 
layered composite as shown in Figure 1 is assumed to contain two 
parallel cracks, one in the mid-plane of the elastic layer and 
the other in one of the half-planes. Figures 2-9 show the varia- 
tions of the four stress intensity factors (at each crack tip) 
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with the material and geometrical parameters entering into the 
problem. First an epoxy layer with elastic properties y-j = 
4.5x10 psi, v-j = 0.35 is sandwiched between two aluminum half- 
planes: y 2 = 10^ psi an d = 0.3. Since the primary objective 
of this work was to study the disturbance problem, the input 
tractions were assumed to be uniform uni-axial stresses with zero 
shear component, i.e., in equation (2) 

P] (x) = q 1 (x) = 0 

P 2 (x) = a Q or (12) 

q 2 ( x ) = a 0 • 

The crack in medium 1 was loaded first and, as expected, we get 
negative stress intensity factors for crack 2. The results are 
shown in Figure 2. Figure 3 shows similar results when only the 
crack in medium 2 was loaded. It is clear that the results need 
to be superimposed if both the cracks are loaded simultaneously. 

In Figures 2 and 3, the stress intensity factors K j, K n» K j 

2 

and K ^ with respect to the distance of crack 2 from the inter- 
face are shown. The two cracks have been assumed to be of equal 
length. We observe that the absolute magnitudes of all the K 
values increase as crack 2 nears the interface. The layer thick- 
ness here was assumed to be equal to the crack length. Figures 
4 and 5 show similar plots for a layer of double the thickness. 
Notice the decrease in the interaction between the two cracks 
since they are farther apart now. Also, the stress intensity 
factors at crack 1 due to the loading at this crack increase due 
to increasing the layer thickness. In limit when h^-~, j = 1.0, 
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12 2 

K ii = K i = K ii = 0. This effect is depicted in Figure 6, 
where the stress intensity factors are plotted with respect to 
the layer thickness. However, a reverse effect is observed at 
the stress intensity factors at crack 2 due to the loading at 
crack 2 itself, as shown in Figure 7. K j and K u decrease as 
the layer thickness is increased and, in the limiting case when 
h-* 00 , these approach asymptotically to the values obtained for 
the problem of a bimaterial medium containing a single crack in 
one of the half-planes. Interaction terms, of course, vanish 
as h-*-°°. Figures 8 and 9 show the effect of the material proper- 
ties of the composite constituents. Here an aluminum layer is 
bonded between two epoxy half-planes. Similar observations are 
made when the crack in medium 2 approaches the interface and 
when the layer thickness is increased. 

As a second example, the case of the elastic layer contain- 
ing two collinear cracks at the mid-plane and located symmetri- 
cally is considered. Since the problem is symmetrical, stress 
intensity factors at only one of the collinear cracks need be 
computed. In all numerical cases, the location of crack 2 and 
the layer thickness have been kept fixed. Again either the 
collinear cracks or crack 2 is loaded at a time. Variation of 
the mode I stress intensity factor with respect to the distance 
between the collinear cracks is shown in Figure 10. When these 
collinear cracks are far away, only a little interaction between 
the parallel crack is observed. Another interesting phenomenon 
observed is that, if the two collinear cracks are close by, i.e.. 
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c/a <1.25 (see insert in Figure 10 for the notation), the un- 
stable crack propagation would first take place at the inner 
crack tip, thus generating a single crack in the layer. On the 
contrary, if the two cracks are far apart, i.e., if c/a> 1.25, 
the outer crack tips would become unstable before the inner 
crack tip. Figure 11 illustrates the mode II stress intensity 

factors for the same problem. When crack 2 is loaded, the 

2 2 

stress intensity factors K j and K jj remain practically unaf- 
fected due to the displacement of the collinear crack locations 
(Figure 12). However, the interaction between these cracks is 
quite strongly affected by the distance between the two coll in- 
ear cracks, as shown in Figure 13. Again, when the two cracks 
come closer, Kj at both the crack tips increases monotonical ly , 
with faster rise in the inner crack tip. Kjj at the inner tip, 
however, undergoes a maximum value and steadily decreases to a 
very low value. 

In conclusion, whenever there is a structure containing 
multiple cracks, an analysis of the type described in this paper 
is essential in order to find the critical configurations under 
which the structure may be most vulnerable. 
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Figure 2. Stress Intensity Factors vs. h -j / a for h/a = 1.0, 
P 2 ( x ) = -a and p^x) = q^x) = q 2 (x) = 0. 






Figure 3. The Case with q 2 ( x ) = an< ^ 
p 1 ( x ) = p 2 (x) = q -j ( x ) = o. 


1.2 



Figure 4. Stress Intensity Factors vs. h -j / a for h/a = 2.0, 
P 2 (x) = -a and p^x) = q^x) = q 2 (x) = 0. 
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The Case with q 2 ( x ) = -a and 
P-j ( x ) = p 2 (x) = q ] (x) = 0. 


Figure 5. 






2.0 - 5 

jjl\ = 4.5 x 10 psi . Z/| = 0.35 
/i 2 = I0 7 psi. 1/ 2 = 0.30 



Figure 7. The Case with q 2 (x) = “ 0 and 
P-, ( x ) = P 2 (x) = q 1 (x) = 0. 



Figure 8. Stress Intensity Factors vs. h-j/a for the Aluminum 
Layer Sandwiched Between Epoxy Half-Planes, 
h/a = 1.0, p 2 ( x ) = -o and q 2 (x) = 0. 



y 


1 

7 

^ 2 =IO pSI. 

V 2 = 0.30 

1 

Mi = 4.5 XIO 5 psi. 

h l- — 2c- 

D Z/| = 0.35 

— 



Figure 10. Stress Intensity Factor Kt vs. the Distance Between 

the Collinear Cracks (a/c). P 2 (x) = -a and q 2 (x) = 0. 



11 K vs. a/c. 
Figure 11 • v c 
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Figure 12. Kj vs. a/c. q 2 ( x ) ~ _a and P2^ x ^ 0 




Figure 13. Kjj vs. 


a/c. q 2 (x) = -o 


and P 2 ( x ) “ 


